Slow light in molecular aggregates nanofilms 
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We study slow light performance of molecular aggregates arranged in nanofilms by means of 
coherent population oscillations (CPO). The molecular cooperative behavior inside the aggregate 
enhances the delay of input signals in the GHz range in comparison with other CPO-based devices. 
Moreover, the problem of residual absorption present in CPO processes, is removed. We also propose 
an optical switch between different delays by exploiting the optical bistability of these aggregates. 
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The optical engineering of the speed of light plays an 
important role in the development of all-optical devices 
for telecommunications. Among the different mecha- 
nisms exploited to obtain slow light P, Q , coherent pop- 
ulation oscillations (CPO) deals with two-level systems 
and is feasible at room temperature 0, |3l- It is worth 
mentioning that there is some controversy in this regard, 
since most of the experiments can also be explained by 
saturable absorption (SA) 0, @] . Being the source spec- 
tral width larger than the coherent hole of the absorption 
profile, some reasonable doubts raise about the existence 
of coherent population oscillations. From a theoretical 
point of view, a rate-equation analysis does not distin- 
guish between both processes, and a density matrix for- 
malism is required. To develop photonic applications, 
a large fractional delay (time delay normalized to the 
pulse length) is desirable in compact devices. Further- 
more, bandwidths up to the GHz or THz-range and a con- 
strained distortion of the input pulses are needed to inte- 
grate slow light mechanisms in nowadays communication 
networks. However, the delay decreases with the signal 
bandwidth, which makes difficult to achieve these objec- 
tives. Moreover, the residual absorption present in CPO 
processes leads to a greatly diminished output intensity 
when longer delays are sought. In this work we propose 
a new optical CPO-based device considering nanofilms of 
linear aggregates of dye molecules, so-called J-aggregates, 
which leads to a fractional delay up to 0.33 for 135-ps- 
long pulses while the ratio between the standard devia- 
tion of the output and input pulses reaches a value of 2. 
More interestingly, the signal modulation is amplified in- 
stead of absorbed, which establishes a clear signature to 
disthinguish between CPO and SA processes in typical 
two- level systems, since in the latter no gain in the weak 
intensity modulation can be achieved. We also show how 
the already predicted optical bistability on J-aggregates 
01 can provide an all-optical switch between two well 
differentiated time delays and pulse distortions. Remark- 
ably J-aggregates display narrow absorption bands red- 
shifted with respect to those of the isolated molecules due 
to collective interaction, which leads to interesting opti- 
cal phenomena studied in the last decades Q . Because of 



disorder effects affecting these systems, only some of the 
molecules of the aggregate are coherently bound. In the 
most favorable situation, namely, at low temperatures, 
the number of molecules over which the excitation is lo- 
cahzed is ~100, although a real aggregate is consisting 
of thousands of individual molecules. In view of these 
properties, we consider J-aggregates as modeled by an 
ensemble of inhomogeneouslv broadened two-level sys- 
tems at low temperatures [3, Q. With regard to CPO 
processes, we study the response of an ultrathin film of 
oriented linear J-aggregates to a strong pump field Eq 
at frequency w, and two sidebands £^±1 at frequencies 
Lj ±5. Here 5 is the beat frequency between fields. Due 
to disorder effects, the film can be considered as con- 
sisting of homogeneous aggregates of different coherent 
sizes N . By using the density- matrix formalism under 
the rotating wave and slowly varying amplitude in time 
approximations, we describe the state of a segment of 
size TV as: 
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Here is the slowly varying in time coherence between 
the two levels which depends on the segment size N. The 
transition frequency and the dipole moment between the 
ground (a) and upper (b) energy levels of every segment 
read and = d^VN respectively. The superscript 
1 refers to single-molecule properties. We assume field 
polarization directed along the transition dipole moments 
of all the aggregates, which in addition are parallel to 
each other as well as to the film plane. The relaxation 
time of the population inversion = (p^ — p^) due to 
spontaneous emission is T/^ = /N, while T2 is related 
to other dephasing processes. 

Field propagation in dense media can be described by 
means of an integral equation where the slowly varying 
approximation in time is considered but not in space 
[9]. For film thickness smaller than the optical wave- 
length, spatially homogeneous polarization can be as- 
sumed, which leads to the following field equation. 
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where /xq is the permeability constant, c is the speed of 
light and L is the film thickness. The incident field is 
= Eq+ Eicyip{-iSt) + E_ieyii>{iSt). The sec- 
ond term accounts for the field created by the molecules 
polarized by the incident field. We checked the full prop- 
agation integral to raise the same results as Eq. for 
films of tens of nanometers, a thickness achievable by the 
spin coating technique Thus, for simplicity and to 
gain a deeper analytical insight, we will restrict ourselves 
to this case. Equation ([2]) resembles that found by Lu et 
al., who proved that local field effects can improve the 
slow light performance in a hybrid nanocrystal complex 

[HI- 

The polarization is calculated by considering the con- 
tributions of all coherent segments of different sizes, 
P = NoJ2NPiN)d^'^^- Here Nq is the density of lo- 
calization segments and p{N) refers to the disorder dis- 
tribution function over localization lengths. Note that 
the size dispersion of the coherent segments in the sys- 
tem results from the inhomogeneous broadening affecting 
the J-band at low temperatures, which mainly gives rise 
to the fluctuation of transition energies . As in Ref. 
01, we replace the average over sizes by one performed 
over the normalized detunings = (cj — aj^)T2. The 
distribution p{N) indeed can be accessible by absorption 
experiments and is considered as Gaussian-like hereafter: 
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Aq = {uj — ujq)T2 being the detuning between the incident 
frequency and the mean of the transition frequency distri- 
bution ujq. The magnitude of the J-band width resulting 
from the inhomogeneous broadening is denoted by G in 
units of 1/T2. From now on size dispersion effects are re- 
stricted to these detuning effects in our calculations, for 
the sake of simplicity. Thus, we will substitute the size- 
dependent quantities by its main value in the aggregate 
and remove the index N, i. e., = Ti and = d. 

Similarly to Ref. [ll], we treat Eqs. ([T]) to all orders in 
the strong field i?o, while keeping only first-order terms 
in the weak fields E±i. Within this approximation, the 
solutions to Eqs. ([!]) are found by considering the Flo- 
quet harmonic expansion: 
a^^expiiSt) and = Zg + Z^exp{-i5t) 
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where ^ = ST2 and r — T2/T1. The DC response of the 



population is Z^c, and Zi accounts for the coherent pop- 
ulation oscillations, which leads to the absorption dip. 

We will refer to the Rabi frequency defined in units of 
I/T2 as f2o.il = 2dEQ ±iT2/h from now on as: 



N 



,±1 



(5) 



Here 7^; = ^Q\d\'^ Nqclu LT2 / 2h is the collective superra- 
diant damping of an ensemble of two-level molecules. 

Once we algebraically solve Eqs. (j?]) and taking into 
account Eqs. ([S|) it can be demonstrated that: 
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where = {^ + ir){^+A'^ +i){C-A^ + i)-\no\H^+i)- 
Equation ^ has one or three roots depending on the pa- 
rameters. This leads to bistable solutions for the strong 
field rio when is larger than a threshold value Q- 

We consider an incident field such that il™ — , 
leading to a sinusoidally modulated intensity with fre- 
quency 6. Thus, Eqs. d?]) and dS]) can be simplified and 
the resulting fields fulfill ili = The transmittance 

T and the dephasing (p induced by the film is calculated 
by the ratio between the output and the input signals: 



- = Te"l', 
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Thus, the fractional delay is defined as F = 0/27r. 

We first analyze the case in absence of size dispersion. 
We take the parameters of Pseudoisocyanine-Br (PIC- 
Br) as it is one of the most studied J-aggregates. Hence, 
we use Ti = 37 ps (corresponding to a homogeneous ag- 
gregate of size TV = 100) and T2 = 0.02Ti. This value is 
consistent with measurements at low temperatures [l3l | 
and allows direct comparison with previous CPO works 
14| . The transition dipole moment is d^ = 12.1 D and 
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the concentration of aggregates is Nq ~ 10^"^ m~'^. These 
parameters give values for 7^ around 10 for film lengths 
of tens of nanometers. As in the usual CPO case, the 
maximum delay is obtained when the pump field reaches 
the saturation intensity, which depends on 7^. Here- 
after, and for each case shown, we will restrict ourselves 
to this optimal intensity. These values correspond to a 
photon flux of ~ 1.4 X 10^^ phot/cm^ for a 100-ps- length 
pulse, and a surface density of monomers ~ 10^'^ cm~^. 
This means that the typical number of absorption and 
re-emission cycles per aggregate is ~ 1, a much smaller 
value than those usually occurring in J-aggregates exper- 
iments [3], and that allows to ensure photostability of 
the samples. Moreover these intensities are low enough 
to neglect the blue-shifted one-to-two exciton transitions, 
which supports the use of a two- level model 0. 

Figure [T] shows the fractional delay and transmittance, 
calculated with Eq. ([9]), as a function of the normal- 
ized modulation frequency and for different values of 7^. 
An increasing value of this parameter, which accounts 
for the collective interaction within the aggregate, gives 
rise to higher fractional delays reaching values of ^0.2 
for 7_R ~ 8. Values of 7_r > 8 allows bistability to ap- 
pear, which introduces a large signal distortion. There- 
fore, only cases up to = 10 are shown. Figure [T] 
(right panel) shows how the transmission increases with 
7i{ as well. It must be noted that the system presents 
gain (T > 1) in the weak probe electric fields for modu- 
lation frequencies close to those revealing the maximum 
fractional delay. We explain this behaviour as a strong 
scattering process from the pump to the weak field due to 
the transient temporal grating generated in CPO. This is 
clearly relevant as CPO processes in other media usually 
present non-desirable residual absorption. Hence, in such 
cases it may be necessary to amplify the signal after prop- 
agation through the slow light medium, specially for long 
media or high ion-densities [l6| . The presence of gain in 
these films could overcome this drawback of CPO-based 
slow light while is revealed as a clear difference with SA 
processes in typical two-level media. 

More relevant to telecommunication applications is the 
performance of the system with input pulse signals. To 
analyze the fractional delay and distortion of the output 
pulses after passing through the ultrathin film we numer- 
ically integrate Eqs. ([T|) together with Eq. ([2]). To vali- 
date the simulations, we first check the numerical integra- 
tion with the results previously obtained. Figure [1] shows 
the perfect agreement with the analytical predictions for 
the case of a sinusoidal input signal. Figure [2] depicts 
the fractional delay {F) and distortion (D) as a func- 
tion of the pulse temporal width (defined as FWHM). 
The distortion is defined as the ratio between the out- 
put and input-pulse standard deviations. Although the 
maximum fractional delay is accompanied by a large dis- 
tortion, values up to -F = 0.33 are obtained with dis- 
tortion 2 for 135-ps-long pulses, which gives ~7.5-GHz- 
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FIG. 1. Fractional delay and transmittance of a sinusoidal 
input signal as a function of the dimensionless detuning ^ = 
&T2 for different values of 7_r and G = 0. The inset shows the 
Output-Input curve for the strong field Vi^. The arrow points 
the optimal intensity for the case 7_r = 6. Solid dots result 
from numerical integration of Eqs. U]) and ([2} . 

bandwidth. These results notably improve previous data 
obtained by CPO-based slo w" lig ht in semiconductor ma- 
terials at GHz-bandwidths [ITTl. As mentioned before, 
increasing values of 7^? lead to larger delays. However, 
the distortion generated when bistability occurs limits in 
practice the values of 77? up to 8 for G = (no size dis- 
persion). Note in Fig. [2] (inset) how the transmittance of 
the electric field is larger than 1 for most of the analyzed 
input-pulse widths, in agreement with results of Fig. [TJ 
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FIG. 2. Fractional delay and distortion for pulsed input sig- 
nals against the initial pulse temporal width. Different values 
of 7h are considered. Inset: Transmission vs FWHM. 

Let us now analyze the influence of the size dispersion 
on the slow light performance. To numerically integrate 
Eqs. ([T]) and ^ with the inclusion of the inhomogeneous 
broadening, we carefully choose the sampling under the 
curve defined by Eq. ([3]) to reproduce the analytical re- 
sults for a sinusoidal signal. We now focus on the ef- 
fects of size dispersion on pulse propagation. As it is 
well known, the presence of inhomogeneous broadening 
reduces the slow light performance, as can be seen in Fig. 
[3l where the fractional delay and distortion are plotted 
against the input-pulse width. A value of G = 3 already 
reduces the fractional delay up to 5 times with respect 
to the value obtained without size dispersion. However, 
as shown at right panels of Fig. |31 the detrimental effect 
of a larger inhomogeneous linewidth can be compensated 
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by increasing the value of 77^. Larger values of can 
be obtained by modifying the temperature or increasing 
the aggregates concentration in the sample Q- 
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FIG. 3. Fractional delay and distortion for input pulses vs the 
initial pulse temporal width (FWHM), for different values of 
G and (a) for -yn = 7, (b) for jr = 7 (G=0), 7h = 12 (G=l), 
■jR = 22 (G=2),7fl = 34 (G=3),7fl = 49 (G=4). 

We finally propose a mechanism to take advantage of 
the optical bistability present in J-aggregates nanofilms 
for slow light applications. Bistability allows to rapidly 
change between two different output intensities for the 
same input. As the fractional delay depends on |rioPj 
this property turns out into a switch between two well- 
differentiated delays and distortions for the same |f^o"P- 
Fig m depicts this process for a value of 7^; = 10 and no 
size dispersion, where the bistability loop can be found 
for 0.69 < i^l'o"f < 0.76. Starting in the lower branch 
of the loop (position (1) in Fig 14]), and after propagat- 
ing the first pulse, a short pulse in the input signal causes 
the system to switch to position (2), in the upper branch. 
Fractional delays achieved in positions (1) and (2), for a 
111-ps-length pulse are 0.43 and 0.09 respectively, while 
the distortion values are 3.6 and 1.07. The lowest switch- 
ing time is limited by Ti, acording to simulations, being 
of some tens of ps for J-aggregates films. 
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FIG. 4. Switch between different delays through the bista- 
bility loop for a pulse FWHM = 111 ps. Arrows mark the 
change of jfiol when applying a short pulse in the input sig- 
nal. In position (1) F = 0.43 while F = 0.09 in (2). Output 
and input pulses in (1) and (2) are shown at the right panel. 



In conclusion we showed that J-aggregates ultrathin 
films can produce large fractional delays by means 
of CPO processes, even in presence of inhomogeneous 
broadening and thanks to the cooperative behavior of 
the aggregate molecules. This system does not suffer 
from residual absorption in the weak probe fields, in op- 
position to the usual CPO-based slow light. We also 
demonstrated how optical bistability could be used to 
produce a fast switch between different delays. We be- 
lieve these organic compounds present a viable alterna- 
tive to semiconductor slow light devices in the nanoscale. 
To this end, our results can motivate further research 
at room temperature operation and telecommunication 
bands. In this sense, the development of new aggregates 
such as Porpho-cyanines [l8| provides great opportuni- 
ties. Moreover, new studies including processes such as 
one-to-two exciton transitions and exciton-exciton anhili- 
tion are currently in progress. These effects seems to 
counteract the killing action of the inhomogeneous broad- 
ening [19l] , which would impose better experimental con- 
ditions to obtain the delays shown in this work. 
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